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Abstract 



^ . The tunneling rate of a skyrmion in ferromagnetic spin- 1/2 Bose condensates through an off- 

, centered potential barrier is calculated exactly with the periodic instanton method. The prefactor 

is shown to depend on the chemical potential of the core atoms, at which level the atom tunnels. 
■ Our results can be readily extended to estimate the lifetime of other topological excitations in the 

o , 

^ O^ , condensate, such as vortices and monopoles. 
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I. INTRODUCTION 



Macroscopic quantum tunneling(MQT), the tunneling of a macroscopic variable of a 
macroscopic system, has recently received much attention in studies of Bose-Einstein con- 
densation(BEC). The tunneling of a condensate through an optical lattice potential |1|, 0] 
provides an atomic physics analogue of a Josephson junction array, while in principle the 
analogue of a single junction can be realized by two condensates confined in a double well 
potential |^. The recent experimental success in all-optical trapping of an atomic conden- 
sate opens the prospect of studies into the internal structure of spinor BECs, including 
the possibility of creating some topological excitations ||^ such as skyrmions, monopoles, 
merons or axis-sjTumetric or non axis-symmetric vortices both for antiferromagnetic and 
ferromagnetic condensates. Among various topological structures, the Mermin-Ho (MH) 
and Anderson- Toulouse (AT) [Q coreless non-singular vortices are demonstrated to be ther- 
modynamically stable in ferromagnetic spinor Bose-Einstein condensates with the hyperfine 
state F = Skyrmions, which do not have an ordinary vortex core due to the spin 



degree of freedom, are also proposed in the spinor BEC |T^, [11], |T2[and are shown not to 
be thermodynamically stable objects without rotation [|13|, |T^. Once created, the radius of 
such a skyrmion shrinks to zero so that one must detect and manipulate it in the duration 
of its lifetime. 

The skyrmion texture in a ferromagnetic spinor condensate can be described conveniently 
by a position-dependent spinor 

C(r) = exp|-l^r.slc^. (1) 



The constant spinor C,^ is the usual basis that diagonalizes the Sz component of the spin 
matrices S and uj{r) is a real function of radius r satisfying the boundary conditions C(j(0) = 
277 and limr_^oo = 0. For the skyrmion with size of the order of the correlation length 



^ = 1/ i/Svranoo or less, where is the average atomic density and a the s-wave interatomic 
scattering length, the problem can be reduced to a nonlinear Schrodinger equation by an 
ansatz for u{r) with the gradient term |VC(r)|^ in the Gross-Pitaevskii energy functional 
regarded as some external potential V{r) = \'VC{r)f /2'm. In the spin-1/2 case, for 
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example, the latter takes the form 113 



. ^ 32 jr/Xf [3 + 2 (r/A)^ + 3 jr/Xf] 
2m [l + (r/A)Y 

for an ansatz uj{r) = 4 cot^^[(r/A)^], where the variational parameter A corresponds phys- 
ically to the size of the skyrmion. The lifetime of the skyrmion is estimated by employing 
a WKB expression for the tunneling rate T = ^e~'^=/^, with Sc the action through the 
barrier and uq the characteristic frequency of the harmonic potential which was used to ap- 
proximate the potential V{r). Due to the inaccuracy of the prefactor in the tunneling rate, 
which makes it difficult to give a reliable result for the decay rate, more efficient methods 
are needed for the investigation of this problem. 

The instanton method as a powerful tool for dealing with quantum tunneling phenom- 
ena has generally been used in the evaluation of the splitting of degenerate ground states 
or the escape rate from metastable ground states ||T5[. A method of evaluating quantum 
mechanical tunneling at excited energy states has been developed recently by means of pe- 
riodic instantons and bounces |16|, which are characterized by nonzero energy and satisfy 
manifestly nonvacuum boundary conditions. Solvable models include, level splittings for the 
double well and sinh-Gordon potentials, decay rates for the inverted double well and cubic 



potentials, and energy band structures of the sine-Gordon and Lame potentials [IG, 17 



The off-centered potential barriers serve as another class of physical systems which permit 
analytical evaluation. 

In this paper, we investigate the tunneling behavior of the skyrmion from the core to 
the outer region through an off-centered barrier. We first solve the equation of motion in 
the Euclidean version to find the classical configuration, which in our case is a bounce. 
In Section III we present the formalism of the periodic instanton theory for tunneling and 
calculate the decay rate exactly. The results obtained are applied to estimate the rate of 
shrinking of the skyrmion in the two-component ferromagnetic Bose-Einstein condensate. 
Finally we summarize the main results. 
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II. INSTANTONS FOR OFF-CENTERED POTENTIAL BARRIER 



In this section we consider the instanton solution for the tunnehng in an off-centered 
potential barrier as depicted in Figure 1 in which the potential 

^^''^ ~ (l + 5r2)2 

takes into account all the main features of the real barriers in the skyrmion excitation in 
the condensate, both in spin-1/2 and spin-1 condensates, for any reasonable ansatz loir). 
We observe here some essential conditions for this simplified model: First, as a function 
the ansatz for lo should decrease monotonically from 27r to 0, since this will correspond 
to the smallest gradient energy for the spin deformations; correspondingly, this excludes 
any oscillation in the decrease of the potential Vir) when r tends to +00. Furthermore, 
Vir) should be an off-centered potential barrier with a maximum height V(rm) at r = r^, 
and 1^(0) = l^(+cx3) = 0. Finally, the potential should be an even function of r, and to 
avoid the point r = becoming a singularity, we have V^"(0) > so that the harmonic 



oscillation frequency ujq can be well defined as A/V""(0)/m. The barrier (^) is just the 
simplest form fulfilling the above requirements, with parameter A determining the barrier 
height, parameter B the position of the barrier: 

Tm = VT/B, = UJo = ^/2A/m. (4) 

To estimate the lifetime of the skyrmion, we calculate the tunneling rate from the core to 
the outer region through a barrier, the core atoms having a chemical potential /Xcore (hereafter 
abbreviated as /i). The first step of the instanton method is the so-called Wick rotation of 
a phase space corresponding to a transformation to imaginary time r = it. After the 
transformation the Lagrangian is replaced by its Euclidean counterpart 

(9 

The classical solution which minimizes the corresponding Euclidean action satisfies the 
equation 

1 f dr,^ ^ 



which can be viewed as the equation of motion for a particle of mass m with energy — E^i in 
a potential —V. For the tunneling process in the condensate, we assume that the skyrmion 




FIG. 1: The off-centered potential and the bounce configuration in two imaginary time periods. 
For the spin-1/2 ^'^ Rb condensate the potential and the radius are in units of h?/2m^'^ and ^, 
respectively, and the parameters are chosen as A = 96 and B = 1, where the size of the skyrmion 
A used is approximately the corelation length ^ corresponding to 20 core atoms. 

has decreased to a size for which the barrier is so high that the overlap between the core 
atoms and the external atoms is exponentially small. The classical turning points on both 



sides of the barrier can be determined by the relation l^(ri 2) = /i as suggested in ref. [|13 



ri,2 



2B 



(7) 



The reason why we can handle a nonlinear problem by means of a linear equation of motion 
is that we discuss the tunneling behavior in the barrier region where the nonlinear interaction 
is negligibly small. Furthermore the condensate at the ground state can be well described by 
a macroscopic wavefunction with unique phase just as in the single particle case. However, 
there are obvious differences between the BEC tunneling system and the usual one-body 
problem, i.e. the nonlinear interaction contributes a finite chemical potential fi, which 
replaces the integration constant on the right hand side of eq.(0). 
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The classical configuration is a bounce which is the solution of eq. (j^) and can be expressed 
in an implicit form 



/(rc) = ^cT. (8) 
Here we have assigned a characteristic frequency 



^ V ^. (9) 



and the function / takes the form 



fi'i^c) = — Ui + Br2 [E{ui) — k'^ snuicdui\ (10) 

where sn, cd are two Jacobian elliptic functions, Ui = F{if, k) and E(ui) are the first and 
second kind of incomplete elliptic integrals with modulus k = — r\lr\ respectively [^], 
and 



^ = sin-^ ^ - (11) 



r2(r2 - rfj 

The solution is subject to the following boundary conditions 



r = 0, r = ri, 

r = ±T,r = ri, (12) 
T = ±T/2,r = r2, 



and exhibits periodic oscillation with imaginary time period 

T = — (-K{k) + Br2E{k)] , (13) 

where K{k) and E{k) are the first and second kind of complete elliptic integrals with mod- 
ulus k, respectively. In Figure 1 we depict the periodic oscillation of this pseudoparticle in 
two periods. A remarkable feature of this bounce configuration is that there is no vacuum 
analogue as in the case of the simple cubic metastable potential, the latter describing the 
tunneling behavior of a particle located at the ground state. As the energy E^ (or the chem- 
ical potential fi) approaches zero, the barrier will become infinitely thick and the particle 
confined in the core region will be stable, with no possibility to tunnel to the outer region. 
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III. EXACT CALCULATION OF THE DECAY RATE 



The tunneling rate of the condensate core atoms was given by a simple expression of the 
form r = Pe~^/^, where P and W are coefficients which depend on the detailed form of the 
metastable potential. The quantity W appearing in the exponential is the Euclidean action 
of the bounce solution and gives the dominant contribution to the tunneling rate, while the 
prefactor P originates from the fluctuation around the classical configuration. For a rather 
rough estimate, P is often taken to be the attempt frequency uj^j^-n as was done in ref. [IB 



However, as we will show below, this simple evaluation is not accurate. This paper provides 
a powerful instanton tool to obtain this prefactor. 

We recall for the sake of convenience the main ideas of the periodic instanton approach. 
Let's first denote the wavefunction of the core atom condensate with chemical potential 



/i by I i)^ >, where | >= ^n{j) originates from the density and satisfies the 
Gross-Pitaevskii equation for single component 

H\%>=fi\%>, H = -^ + Vir)+g\^,\' (14) 

The effective external potential V{r) [|14| comes from the gradient term of the spinor | VC(r) f, 
and the term with coupling constant g represents the strength of the interatomic interactions. 

The tunneling effect leads to the decay of the metastable state. In the case under discus- 
sion, the nonconservation of an exponentially small probability current through the barrier 
requires that the chemical potential has an imaginary part proportional to the decay rate 
P^ , r = |lm/i. Consider the transition amplitude from the state \ ip^ > to itself due to 
quantum tunneling in Euclidean time period T. The amplitude is simply 

A =< il)^\e-"^'^ I >= e-^^l''. (15) 

In general the transition amplitude is calculated with the help of the path integral method 

as 

^ = y" '^*^irf)'tp^{ri)K:{rf, T; r^, 0)dridrf, (16) 

where Vf = rc{T),ri = rc(0) denote the end points of the bounce motion, which tend to the 
turning points ri (see the boundary condition eq. (p!2D). The wave functions 'ip^{ri),ip^{r f) 
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in the barrier region are specified in the WKB approximation as pi 

C 



P 



exp 



v/2m(/i - 1/(r)) 



pdr 



(17) 
(18) 



with C a normahzation constant to be determined below. The Feynman kernel is defined as 
the summation over all possible classical paths r(r) 



/C(r;,T;r„0) 



V{r}exp{-S/h). 



(19) 



We know that the classical solution (|^) which minimizes the action S gives rise to the 
major contribution to the above kernel integral, while the quantum fiuctuation around it 
results in a prefactor P. In the period T the bounce eq. (||) completes one oscillation and 
crosses the barrier region twice, back and forth. The Euclidean action is thus calculated in 
this period as 



/■^ / fdr. 
C{r, r)dT = \ m 



dr 



+ yU 1 dr 



(20) 



while the so-called abbreviated Euclidean action pO| 



W 



r2 

dr^/2m(y{r) — jj,) 



ri 



can be expressed in terms of elliptic integrals 
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W = — 

UJc 



— ^— — 4( - —K{k) - i2Br2E{k) 



_r2{l + Brf) r2 
with n(Q;^, k) the complete elliptic integral with the parameter 

o Jv 



a 



(21) 



(22) 



(23) 



1 + Brl ■ 

It is obvious from the potential that A, B > 0, so < < k"^ , the third elliptic integral is 
complete and belongs to the case III[|l9| . 

The imaginary part of the chemical potential can be derived by considering the amplitude 
A as the sum of contributions from any number of bounces [|16|. The zero bounce contribution 
results in the real part of the chemical potential 



^(0) 



(24) 
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FIG. 2: The decay rate as a function of the chemical potential. Inset: the chemical potential 
dependent frequency lo. The solid curves represent our exact result while the dotted curves corre- 
spond to the case for constant attempt frequency wq- All curves are calculated for the parameter 
of a Hh spinor condensate and the chemical potential ii is given in units of the barrier height Ym- 

The one bounce contribution comes from the classical configuration with period T, and can 
be obtained by expanding the kernel (|1^) around the bounce (H) 



m 



(25) 



Generalizing to the case of n bounces straightforwardly, i.e., assuming the pseudoparticle 
completing n oscillations in the period T, one has 

nl \m 

The total transition amplitude is given by the sum over all bounce contributions 



g -nW/h^-fiT/h_ 



(26) 



(27) 



The imaginary part of the chemical potential is obtained by comparing eg. ([27|) with eq.(|T5| 



Im /i 



-W/h 



m 



(28) 
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which results in the decay rate 



1 2 



Im /i 



-W/h 



(29) 



2h ' m 

where the factor 1/2 comes from the analytical continuation. Physically, this results from 
the assumption in the decay problem (and not in the MQC problem) that the wave that has 
tunneled will never return. Mathematically, it is due to the fact that the deformed contour 



runs from to ioo, and not from —ioo to ioo[15 



The constant C can be determined from the normalization of the wave function in the 



classically accessible region, which is connected with those in the barrier region through pi 
C 



exp 



— / pdr 



2C 



cos 



(30) 



'\P\ ^ \ ^Jr/ J 

We restrict the integration in the classically accessible region, i.e., in the potential well, 
r < Ti, since outside of this range ip decreases exponentially. Because the argument of the 
cosine in the wave function is a rapidly varying function, we can, with sufficient accuracy, 
replace the squared cosine by its mean value 1/2. This gives 



Inserting this into the decay rate we have 



i = — exp 
27r ^ 



W_ 



where uj is the frequency of the classical periodic motion 

27r TT 



o;(/i) 



2m/ f v^/- 



dr 



and can be calculated as 



IT 

UJr — 



— + 5r2 K(k') - Br^Eik') 
r2 



(31) 



(32) 



(33) 



(34) 



with the complementary modulus k' = Vl — k'^. It must be recalled that the frequency u is in 
general different for different levels, being a function of the chemical potential. We find that 
our expression for the decay rate eq. (|3^) is more accurate than that of Refs. |l^, 14], i.e., in the 
prefactor a chemical potential dependent frequency replaces the constant attempt frequency 
uq. In Figure 2 we show the dependence of this frequency on the chemical potential; it 
decreases from c^o as the chemical potential increases from 0. This factor suppresses the 
tunneling rate greatly when the chemical potential approaches the barrier top as shown in 
the figure, which would be expected to increase the lifetime of the skyrmion. 
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FIG. 3: Ansatze for and the corresponding potential barriers V{p)- The potential and the 
radius are again in units of It? /2m^'^ and ^, respectively. 



IV. NUMERICAL RESULTS FOR SKYRMIONS 



The skyrmions in a ferromagnetic condensate are energetically unstable as shown in ref. 
1^ . The time scale on which the skyrmion shrinks may be evaluated for two cases: For a 



large skyrmion, its size decreases at a rate T large ~ 18 sec~^ ^/A for ^'^ Rh spin-1/2 condensate 
of central density lO^^cm"^ and realistic experimental conditions. For skyrmions with sizes 
of the order or less than the correlation length S,-, the shrinking rate is determined by the 
tunneling rate from the core of the skyrmion to the outer region. In previous studies []T3|, |T^ 
the authors estimated roughly the lifetime of this small skyrmion due to the tunneling process 
employing a WKB expression for the tunneling rate. In this section we thus reconsider the 
lifetime using the result derived above, with the modification originating from the prefactor 
included. 

With an ansatz for uj{p) the problem is simplified to a nonlinear Schrodinger equation 
with the external potential of the off-centered form. For a different functional behavior of 
co'(p) it turns out that the effective potential will not be very different, as long as the ansatz 
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satisfies the boundary conditions 0^(0) = 27r and linip^oo ^(p) = and falls off monotonically. 
Here we compute the decay rates of different skyrmion textures with the same size, by taking 
into account two ansatze which were proposed in Ref. |T^ as trial functions for simplifying 
the pair of nonlinear and coupled equations (eq. (8) and (9) in ref. |T^). Namely 

a;i(p) = 4cot-i[(p/A)2], (35) 
^2(P) = T^^' (36) 

where p = rj^ and A corresponds to the size of the skyrmion and is also given in units of ^. 
Considering the large distance behavior of the coupling equations for n(p) and cc'(p), we see 
that for large skyrmions, the density fluctuations scale as and so should the ansatze 

for ioi^p). This is the reason why we would arrive at a non-physical result for a seemingly 
reasonable ansatz ujj,{^p) = 27rsech(p/A). We also check that non-monotonic behavior, i.e.. 



oscillations in the falling of uj{p), for example, 0^4 (p) = 27r 



will inevitably lead to 



p/A 

singularities in the density profile, though the effective potential holds an off-centered form. 
We show in Figure 3 these ansatze for uj{p) and their corresponding effective potentials V{p). 

One important parameter we should determine is the chemical potential of the core 
atoms because we should know at which level the atom will tunnel out. In principle one 
should solve the two-coupled nonlinear differential equations and derive the density profile 
and spinor(or the function uj{p)). As mentioned already above, we employ alternatively a 
simple approach, i.e., by introducing the ansatz for uj{p). Then from the resulting density 
distribution we calculate the energy for a particular value of A, then the core chemical 
potential can be calculated numerically by differentiating the energy with respect to the 
number of core atoms. Performing the calculation within a Thomas-Fermi approximation, 
which means in the expression for the energy we neglect the kinetic energy term , we finally 
obtain the chemical potential of the core atom p for different values of A and ansatz. 

For the ansatze ui and uj2 we calculate the corresponding chemical potential for A = ^. 
The shrinking rates of the corresponding skyrmions are calculated according to our decay 



rate expression eq. (|32D with the action given by eq.(pi|) and the prefactor given by eq.(p^). 
Figure 4 gives the tunneling rates as a function of the number of core atoms. The calculation 
was performed for a ^"^Rb spin-1/2 condensate with a scattering length of a = bAnm. 

We observe that [22] the correction resulting from the accurate prefactor uj{p) for uoi is 
minor but significant for uj2- In Figure 2 we could generally take the range of the chem- 
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FIG. 4: The shrinking rate of skyrmions as a function of the number of core atoms for uji and 
u}2- The calculation was performed for a Rb spin-1/2 condensate with a scattering length of 
a = 5Anm. The dashed lines are the WKB calculations in Ref. [jl3|, while the solid lines show our 
results from the periodic instanton method. 

ical potential from to Vm- Unlike the situation in a harmonic trap where the chemical 
potential /i increases with the number of condensed atoms as N"^/^ in the Thomas-Fermi 
approximation, in our case fi decreases with N instead. This is because of the fact that the 
trap frequency (for uJi) ujq is inversely proportional to the equilibrium skyrmion width Aq, 
which in turn increases with N(apparently faster than N^^^). Here Aq is determined from 
minimizing the total energy taking into account the outer region of the skyrmion. This 
restricts us to a special domain of /i. For numbers of core atoms ranging from 1 to 20, fJ^/Vm 
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ranges roughly from 0.16 to 0.04. In this interval uj{ij)/ujq starts from 0.985 for Ncore — 1 
and ends at 0.998 for Ncore — 20, which are almost indistinguishable in Figure 4. However, 
for ansatz uj2 the chemical potential (in units of /2m^'^) starts at 42 for one core atom 
and ends at 9.3 for 20 core atoms. The corresponding correction is shown in Figure 4. 

Till now there is still no clear experimental evidence for the skyrmions in the condensate. 
From the above calculations we see that the result for the decay rate depends crucially on the 
detailed form of the ansatz uj. It remains a challenging task to solve the coupled nonlinear 
equations numerically, and to compare the results with those above. 

V. CONCLUSION 

We present here an accurate calculation of the tunneling rate for a class of off-centered 
potentials with a periodic instanton method. Apart from its application to the study of 
the stability of the skyrmion excitation in the two-component ferromagnetic condensate, the 
bounce for the off-centered potential barrier is itself a novel configuration from the viewpoint 
of the scalar field theory. The exact prefactor of the decay rate has been calculated and 
we found it depends on the chemical potential at the level of the atoms tunneling to the 
outer region. This modifies the result for the rough estimate of the lifetime by a constant 
attempt frequency Uq. One can easily find some similar off-centered potentials in other 
topological excitations such as vortices, monopoles, etc. Our periodic instanton formalism 
can be extended to the investigation of the lifetime and tunneling behavior in these systems. 
Further studies should include the properties of the quantum-classical transition of the decay 
rate when the chemical potential increases and surpasses the barrier height. 
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